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ABSTRACT 



We present a geometrodynamical method for determining distances to orbital streams 
of HI gas in the Galaxy. The method makes use of our offset from the Galactic centre and 
assumes that the gas comprising the stream nearly follows a planar orbit about the Galac- 
tic centre. We apply this technique to the Magellanic Stream and determine the distances to 
all points along it; a consistency check shows that the angular momentum is approximately 
constant. Applying this technique to the Large Magellanic Cloud itself gives an independent 
distance which agrees within its accuracy of around 10%. Relaxing the demand for exact 
conservation of energy and angular momentum at all points along the stream allows for an 
increase in orbital period between the lagging end and the front end led by the Magellanic 
Clouds. Similar methods are applicable to other long streams of high-velocity clouds, pro- 
vided they also nearly follow planar orbits; these would allow otherwise unknown distances 
to be determined. 
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1 INTRODUCTION 

High-velocity clouds (HVCs) are defined as HI gas which have 
anomalous velocities that are incompatible with what is expected 
fr om a simple Galactic rotation model. Ever since their discovery 
bv lMulIer et al.l il963l) , their role in the formation and evolution of 
our Galaxy has been a source of speculation. Amongst the attempts 
to explain their origin, the two most well-known and contested are 
the Galactic Fountain model, in which gas is either being thrown 
out of the Galactic disk from supernovae or being heated and ion- 
ized by massive stars, and the alternative model in which the clouds 
are pockets of neutral hydrogen left over from the formation epoch 
of the Galaxy, and which are only now being accreted. 

Given the vast sky-coverage of HVCs, the two theories are not 
necessarily contradictory, as it is likely that some clouds belong 
to each model. However, with the exception of those HVCs that 
are obvious comp onents of a stream with a known progenitor, such 
as the Magellanic i Mathewson et alJI 19741) and Sagittarius streams 
jPutman et al.l l2004'). the lack of distance measurements has hin- 
dered the progress in understanding the role that HVCs play in the 
Galaxy. 

Since newly accreted gas would have a lower metallicity 
than recycled gas from the Galaxy, metallicity measurements, 
if available, should be able to distinguish the ap plicability of 
each theory on a cloud-by-cloud basis (e.g. ISembach et al...200Z ; 
Ivan Woerden et al.ll2000l ; IwakkeilllOol) . but it would still be in- 
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structive to know exactly where these clouds lie. Observationally, 
the difficulty in determining direct distances arises from the fact 
that most of the HVCs do not host stars. If, however, there is one 
star behind and another to the front of the cloud along the line of 
sight, both of which have known distances, then a distance bracket 
can be set for the HVC. The high resolution spectrum of a star di- 
rectly behind the HVC will exhibit absorption lines at the velocity 
of that cloud; features that will be absent in the spectrum of a star 
to the front of the cloud. Starting with the first such measurement 
for Complex M in the early 90's JDanlv etalJl 19931) . other HVC 
complexes f or which a di stance bracket that now exists includ e 
Complex A ( IWakkeil l200l') and Complex WB JThom e t al."2006'). 
The high latitudes of much of the neutral hydrogen gas in question 
means there is a lack of background stars bright enough for high 
resolution spectroscopy, so this technique has not been without dif- 
ficulty. 

Many detailed sim ulations of streams such as the Magel- 
lanic Stream exist (e.g. Connors et al. 200f); Mastropietro et 
2005; Gardiner & Noguchi 1996; Fujimoto & Murai 198 
Fuii moto & Sofud 119761) , however, all such simulations have 
many parameters which must be tweaked to give the best fits to 
the data. In this paper, we describe an alternative set of methods 
with the aim of probing the elusive distances to streams of HI 
gas comprising the HVCs. These methods instead make strong 
hypotheses that enable us to attempt to read the distances out of the 
data themselves; the studies presented here are neither simulations 
nor attempts at modelling. The following section introduces the 
geometrodynamical methods for determining distances to streams 
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Galactic Plane 




Stream Plane 



Figure 1. An orbital stream lying in a plane through the Galactic centre. 
The unit vector cr lies along the stream at a general point P, which is seen 
from the Sun in the direction 1 at a distance d. Its observed radial velocity, 
corrected to the Galactic Standard of Rest, is »;; = v.J; at the position la- 
belled Z, Vi = 0. The unit vector s is along the projection of & onto the sky, 
so s.l = and s is o- — (ct.I)! normalised to unity. At Z, s and & coincide. 



of HVCs and discusses which streams are appropriate for each 
method. We developed the methods in response to a challenge 
laid down by Majewski, who asked whether it was possible to 
determine distances from ra dial velocities al one, when those 
were known along a stream JMaiewski' 2004). In Section 12.21 
we consider distant streams and use heliocentric radial velocities 
corrected to the Galactic Standard of Rest (GSR) as a lowest order 
approximation to the Galactocentric radial velocities. Section |3] 
shows how these radial velocities can be corrected iteratively, 
and shows that the corrections are negligible for the Magellanic 
Stream. The choice of Galactic potential is briefly discussed 
in Section |4] On the assumption that the Magellanic Stream is 
composed of material that was torn off the Magellanic Clouds 
at an earlier passage through perigalacticon, it seems reasonable 
to assume an angular momentum gradient to be present along 
the length of the stream; this topic is addressed in Section |5] We 
present our summary and conclusions in Section|6] 

There is a great body of literature on streams and associations 
of bodi es in plane s through the Galactic centre. The Sagittarius 
stream jlbata et alj|200ll : IPalma et al.1120031) features prominently 
in this literature, however it is not relevant to the methods devel- 
oped here which work off the parallax due to the Sun's displace- 
ment from the plane of the orbit. This displacement is almost zero 
for the Sagittarius stream. 

Some frequently used symbols and their definitions are sum- 
marised in Table[T]with a schematic guide to a generic point on the 
stream given by Figure [T] 



2 DISTANCE DETERMINATION METHODS 

In the following sections, we explore several different variations of 
our geometrodynamical distance-determination scheme and their 
applicability. 



2.1 Near Method 

Let us suppose that we have a set of observed heliocentric radial ve- 
locities relative to the Local Standard of Rest (LSR) along a neutral 
hydrogen stream, which can be re-expressed in the GSR. We define 
a location Z to be a point on any stream where the line-of-sight ve- 
locity as seen in the GSR, vi , is zero. The assumptions made are that 
the clouds of neutral hydrogen form a single stream, the streaming 
material is in a plane containing the Galactic centre, and that the 



stream is moving along itself. If these assumptions are valid, then 
the implication is that the line-of-sight vector to Z lies perpendic- 
ular to the true direction of motion of the stream there. Here, the 
apparent direction will equal the true direction of the stream's mo- 
tion. Hence for a given Z point at any assumed distance of dz 
from the Sun, we will have two vectors that lie in the plane of the 
stream: the vector from the Galactic centre to Z and the directional 
vector, sz, in the plane of the sky. If dz can be found, the plane of 
the entire orbit could be found from the cross-product of these two 
vectors. All other points along the stream could then be projected 
along the lines of sight onto this plane, enabling their distances and 
all three components of their velocities to be found. 

The distance to the Z point may be found by looking for a 
value which fulfills the requirements of energy and angular mo- 
mentum conservation along the stream, assuming that such a value 
exists. 

Let us consider a point on the stream, to which we have a 
line-of-sight vector \j. We denote the unit vector describing the 
apparent direction of the stream on the sky as s. Here and hereafter, 
hats denote unit vectors. If ct is defined to be the true direction of 
the stream at this location and n the normal to the stream plane, 
then by construction, l.s — and a — sees 6* + 1 sin 61, where 9 
is the angle between the apparent direction and the true direction. 
Clearly, (T lies in the stream's orbital plane, and hence 6 satisfies 



cote = 



l.n 



(1) 



The velocity of the stream is, by hypothesis, parallel to a, and 
hence may be written as v = va. The line-of-sight velocity in 
the GSR is then given by vi — v.l = t;sin6. This gives for the 
velocity: 



■ VI 



fi X (1 X s) 



for 



L<t/0. 



(2) 



At the Z point, both vi and 9 are zero. Let us denote quantities 
associated with this point with a subscript Z and let Ro be the 
vector from the Galactic centre to the Sun. If r is the vector joining 
the Galactic centre to a general point in the stream at distance d 
from the Sun, then r = Ro + dl and, specifically at Z, vz ~ 
Ro + dz^z- By definition, r lies in the stream's orbital plane, and 
hence: 

^-" (3) 



= r.fi = Ro.n + dl.n 



d- 



l.n 



and we can eliminate d to obtain an expression for r in terms of 
only Ro, 1 and n: 



r = Ro 



^- _ Ro(Ln) - l(Ro.n) _ (1 x Ro) x n 



Ln l.n 

Equipped with these vector relations, we now show how con- 
serving the angular momentum along the stream enables us to find 
the distance to the Z point, which serves as a representative dis- 
tance to a stream. 

The specific angular momentum (henceforth referred to as the 
'angular momentum') associated with the orbit of a stream is given 
by 



^ The vectors, unless specified otherwise, are stated in the Galactic coor- 
dinates, whose origin is defined as the location of the Sun, and whose x 
direction lies towards the Galactic centre. The z axis points to the North 
Galactic Pole. 
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Table 1. Definition of commonly used symbols, where hats denote unit vectors. 



Symbol Equation where first used 



Description 



Apparent and true direction vectors for a particular point on the stream 

Vector normal to the stream plane 

Unit vectors to a stream point, from the Sun and the Galactic centre 

Heliocentric line-of-sight velocity corrected to the Galactic System of Rest 

Vector joining Galactic centre to Sun = (~Ro, 0, 0) 

Distance to a specific point on the stream from the Sun and from the Galactic centre 

Specific angular momentum and energy associated with the stream 

Ratio of distances to the Z point and the Galactic centre from the Sun 

Galactic potential 

Galactic circular velocity at Ro (unless specified otherwise) 
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Figure 2. A schematic diagram showing the angle ^^ subtended at the 
stream by Ro projected onto n, the normal to the stream's orbital plane. 
Ro joins the Sun to the Galactic centre, i^^j is given by arcsin(Ro.n/ci) 
and is useful in quantifying the range of validity of the near method - see 
main text for details. 



= [(1 X Ro) X n] X [n X (1 X s)] R^, vi / [(l.n)(s.n)] 

= [Ro.(ix s)/(s.n)]i?o«in. (5) 

Conservation of angular momentum implies that the coefficient of 
n is constant along the stream. Since Ro and vi are known, the only 
remaining unknown is n, assuming that s is known at each point of 
the stream from looking at the sky. Since the normal to the plane is 
parallel to rz x sz = (Ro + dzlz) x sz, it may be written: 



n = (rz X sz)/\rz x sz\ = Ro 



[Ro X sz + Dlz X sz] 



rz X Sz 



(6) 



where D — dz/Ro is unknown. On the basis that h is constant, we 
can derive an expression to determine D through other observable 
or calculable quantities. Since h is parallel to n, h.fi is a constant. 
Taking the dot product of h in equation l|5j with n gives an expres- 
sion in which s.n can then be expanded via n given in equation 
l|6j. This resulting expression is then inverted to give the following 
condition along a planar stream: 



[Ro X Sz + Dlz X sz].s 

Ro.(l x s)RoVi 
Setting 



= constant 



Y 



(Ro X Sz) .§ 



and X — — 



(Iz X Sz) . 



(7) 



(8) 



Ro.(l X s)i?oWi Ro.(l X s)7?o -y! 

plotting Y against X should give a straight line of slope — D. 

The offset of the Sun from the Galactic centre produces in- 
sufficient parallax for this method to work if the stream plane 
passes close to Ro. We can see this again from equation ([8](, as 



both the numerators and denominators vanish when Ro lies in the 
stream plane. We performed tests on simulated orbits of varying 
orientations, which were given some random scatter in direction 
cosines in order to test the accuracy to which s was required. We 
found that the angle i?d subtended at the stream by Ro projected 
onto n needed to be at least 9° if errors in the determination of s 
and Sz are not to dominate in equation Q. Figure |2] shows how 
§d = arcsin(Ro.ii/d) is subtended at the stream. It is also true 
that if D is very large, then errors in X may dominate in deter- 
mining the gradient of the line Y + D X — const. Only lower 
limits to D could then be determined. In general, an 'ideal' data set 
will give the expected straight line, with singular behaviour in the 
form of another line (with positive or negative gradient) for points 
which are close to Z. This enables dz to be determined reliably to 
an accuracy of within a few per cent. 

In practice, the data for the Magellanic Stream did not produce 
a straight line. With hindsight we might have expected this, since 
the method should not work for a distant stream for which we have 
a very radial view and very little direct information of the transverse 
velocity components. The radial velocities do not contain enough 
information on the angular momentum of the orbit if the stream is 
distant and, for most parts of the stream, 1?^ subtended at the stream 
becomes too small for s to be determined reliably. Ironically, this 
angle is largest close to Z, where the denominator in equation ^ 
vanishes! The method should be better suited to closer streams but 
there the assumption of a planar orbit is more doubtful unless the 
orbit passes close to the Galactic Pole. We note in passing that the 
plane of the Sagittarius Stream lies much too close to Ro for it to 
be a suitable candidate. 

Our first attempt to use this method on Complex A produced 
rather perplexing results. We found a good straight line between 
the X and Y variables of equation ([Sj, but with a positive slope, 
rather than the negative slope that we would expect. The solution 
gives a distance of — 5 kpc, in other words in the opposite direction 
to which we see it. Clearly this was nonsense. Having tested this 
near method on simulated orbits, we realised there were two possi- 
ble explanations as to why the method might not work for Complex 
A, even if the stream was following an orbit. The angles -ffd at all 
points along the stream are large enough at distances that we might 
expect to find for Complex A. However, there is a competing ef- 
fect that these points are all rather close to the Z point, and the 
denominators in equation ^ vanish at Z. On the other hand, it is 
possible that the stream is simply too close to be in a well-defined 
plane. The flattening of the potential due to the disk would be a 
significant deviation from the assumption of a spherical potential. 

We recently posed the que stion of whether the stellar Orphan 
Stream JBelokurov et alJuOOTD and Complex A might be related 
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jjin & Lvnden-Beli2007l) . By applying the near method presented 
here to the orbit for Complex A found in this work, we found that 
the perplexing positive slope could be explained, in that case, as 
part of the singular behaviour that is expected with points close to 
Z. We also tried the near method on Complex C, which has two Z 
points; see Section |23] for an interpretation of the results. 

In summary, the range of validity of this first method of dis- 
tance determination is to systems with orbital planes well away 
from the Sun-Galactic centre line but for which the angle -dd = 
arcsin(Ro.n/d) subtended at the stream points is at least 9°. This 
leads to distances not more than a factor of order six larger than Ro 
for a stream of Magellanic Stream-like orientation. In Section lAl 
we give a working example of the near method. 

The advantage in attempting to use conservation of just the 
component of angular momentum about the Galactic axis is that it 
requires no detailed model of the Galactic potential in which the 
stream flows. We have, however, identified some conditions under 
which the near method via angular momentum conservation does 
not work. We considered next the implications of energy conser- 
vation, which should not have suffered from the same issues. On 
applying this method to Complex A, we found that the s vectors 
were too inaccurately determined to find a planar orbit which gave 
a consistent distance to the Z point, partly because Complex A ex- 
tends only 50° across the sky with the Z point to one end. Points 
near Z give inaccurate results as they suffer from small denomina- 
tors. 



2.2 Far Method 

Application of the method described in Section 12.11 to the Mag- 
ellanic Stream caused some problems. We need to determine the 
plane of a stream to such an accuracy that the line of sight to the 
clouds can be projected to hit that plane at their true distances. 
Lines of sight at small angles to the plane lead to inaccurate dis- 
tances due to small offsets of the clouds from the best-fitting plane. 
The next method circumvents this by using the line-of-sight vectors 
to individual points on the stream. 

When the points of a stream are at distances from the Sun 
much greater than Ro, then to first approximation the heliocen- 
tric line-of-sight radial velocities corrected to the GSR, vi, will be 
close to r, the Galactocentrito radial velocities. On the hypothesis 
that the clouds of the stream follow an orbit in a plane through the 
Galactic centre, the energy equation in a potential 4>{v) is given by 



-ip{v) =e, 



(9) 



where e is the specific energy (hereafter referred to as the 'energy'). 
If there is a cloud at which the stream reaches a maximum value of 
\vi\ ~ |f|, then labelling such a point the H point and differenti- 
ating equation 1(9}, we find h^ /r'u = —dip/dr\ru = v'i{rH)/rH, 
giving h = rHVc{rH)- This relationship is particularly simple for 
the special case when Vc = constant, in which case a simple log- 
arithmic potential T/'(r) — —v'^ In r/r^ can be used, for some halo 
cut-off radius Vh', we will therefore work with that special case first 
and return to the more general case in Section [4] Subtracting the 
energy equation at a general point from that at H, we have, using 
the above value for h and dividing by v'^/2: 



^ The Galactocentric coordinate system is defined by translating the Galac- 
tic coordinate system along its x axis by Ro so that the coordinate origin 
coincides with the Galactic centre. 



•2-2 2 2 

2 1^ 1 = ^- - 111 ^- ' 



(10) 



where vi^h and vi are used for th and r respectively and Vc is set to 
220kms~^. We may then solve this equation for r/rn- Note that 
there are two possible roots to the equation, one on each side of rn, 
but in practice r < r// in the Magellanic Stream, which eliminates 
one solution. This would clearly change if clouds existed past the 
H point. 

Equation l llOt gives ratios of distances from the Galactic cen- 
tre to points on the stream. The distances from us, d, are then given 
in terms of rn by 

2 

(11) 



(f - 2di.Ro +rI^ ( — 



2 



so for three points Vi on the stream, i=l, 2 and 3, 



di — li.R-o 



Th 



2 .1/2 

^) ^ R^o + {l.-Rof ] . (12) 



The condition that three points with distance vectors d; = diU 
lie in a plane through the Galactic centre is that for these points, 
di + Ro lies in the plane of the stream. Hence we have 



[(di + Ro) X (d2 + Ro)] .(da + Ro) = . 



(13) 



Rearranging this, we can then solve for the value of vh/Ro that 
satisfies 



Ro. (di X d2 + d2 X ds + da x di) 



-1. 



(14) 



(di X d2).d3 

Using equation l llOt . distances may be found to the rest of the 
stream clouds, for which velocity data are available. The mean 
line through the observed velocities as a function of angle along 
the primary stream direction reaches 205 km s^^ some 74° be- 
hind the Z point, which we take to be the H point since it has 
the highest velocity of any substantial hydrogen on that mean line. 
For Ro — 8.5 kpc, a typical set of values found is vh = 75 kpc 
and n = (0.995, 0.080, -0.066). The choice of the three fiducial 
points is somewhat arbitrary, but to get good angular leverage, we 
take two points near the ends of the stream and one toward the 
centre. Our results for the Magellanic Stream are given in Figures 
[3] and |4] and Table [2] We find that the perigalacticon is located at 
approximately 45 kpc and our last point on the stream (just prior 
to the H point) is 70 kpc from the Galactic centre. The centre of 
mass of the Magellanic Clouds is independently determined to be at 
47 ± 5 kpc. The error in determining the distances is approximated 
by the ratio of the uncertainty in the central line of the stream to the 
parallax due to the Sun's offset from the Galactic centre, which is 
about 9%. 

2.3 Two Z Points 

If a stream exhibits two Z po ints, as suggested by the ve locity field 
of Complex C presented by ISchwarz & de Boeii ( l2004h . we have 
yet another possible variation on a theme. 

2.3.1 Method 

Let us denote the unit directional vector for a Zi point as Si. Then, 
by using the fact that s^ for both Z points lie in the plane of the 
stream, we can find the normal n to the plane as 



Si X S2 

§1 X §2 



(15) 
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Table 2. Galactocentric distances to the Magellanic Stream. 
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Figure 3. The Magellanic Stream depicted in its own plane, with the leading 
end (position of the Magellanic Clouds) to the left of the plot. The Galactic 
centre is denoted by the star. The yi and zl axes shown here are nearly 
parallel to the Galactic y and z axes respectively. There is little scatter in 
the third (unshown) dimension, the axis for which lies along the normal to 
the plane. 
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Figure 4. The magnitude of the ^/-component of the angular momentum 
for points along the Magellanic Stream is plotted against distance along the 
stream in the yl direction. An explanation of the error bars is given in the 
main text. 
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Table 3. Columns 1 and 2 give the Galactic longitude and latitude respec- 
tively; column 3 gives the heliocentric line-of-sight velocities connected to 
the Galactic Standard of Rest; column 4 gives the uncorrected distance r 
from the Galactic centre in the plane (i.e. ^Jyl^ + zl^) and the correction 
Sr (see Section [3); columns 5, 6 and 7 together give the position vector in 
the Galactocentric coordinate system. The data coincide with the order of 
points in Figure [5] starting from the left. Note that Sr has been rounded to 
the nearest decimal. The fiducial points used to define the stream plane are 
denoted with f. 



As before, the absolute distance of the plane from us is found by 
constraining the Galactic centre to be in the plane. We know that 
at the Z points, l^.Si = 0, and equation l[3j applies to each of the 
points. The probl em then reduces down to fi nding the s^ vectors. 
From page 31 of Ivan Woerden et al.l 1 120041) . we have an all-sky 
Aitoff projection map of the high- velocity clouds. By eye, it is pos- 
sible to draw what would appear to be the directional line of the 
stream, especially as the stream is so long. We can hence draw a 
directional 'vector' from each of the Z points. Consider the first Z 
point given by Galactic coordinates {Izi, bzi) and a close-by point 
along this 'vector', whose coordinates are given by {Ipi, &pi); si 
will then be given by the vector joining these two points. 

Consider the celestial sphere of radius r = 1. Then si — 
Rlpi — rli where 7? 7^ 1. For the moment, R is an unknown. 
However, since li.si = 0, then by substitution and rearrangement, 
we find that R — (li.lpi)^^. This implies that 



Ipi t 
Si = li . 

Ipi.li 



(16) 



By evaluating the cross product Si x S2 and normalising, the unit 
normal to the plane may be found. 

We now have a way of determining the distance to any point on 
the stream by using equation l|3), which is general and not confined 
to the Z points. It does not assume that the path of the stream is 
confined to an orbit of constant energy and angular momentum. 
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Figure 5. Complex C depicted in its own plane. The Galactic centre corre- 
sponds to the coordinate location (0,0) and the position of the Sun is pro- 
jected onto the plane and shown by 0. The xl and yl axes shown here are 
nearly parallel to the Galactic x and y axes respectively. 



however we require the stream to be in a planar configuration as in 
all previous methods. 



2.3.2 Results 

The velocity field of C omplex C is given by Figure 4 of 
ISchwarz & de Boen J2004h . from which we interpolate the follow- 
ing two Z points: (Zzi,fezi) = (37°, 18°) and {Iz2,bz2) = 
(110°, 48°). With the following choice for associated points to 
find the apparent direction vectors for each of the Z points — 
(Zpi,fepi) = (13°,11°) and (;P2,&P2) = (90°, 42°)— we 
find the normal to the supposed plane of this complex to be 
(0.252, —0.071, 0.965) and a set of distances to the stream with 
a range of '^ 3 — 6 kpc under this method, as shown in Figure 
[S] The stream has a significant width, across which there is some 
non-negligible scatter in the velocities. 

It is slightly concerning that some of the clouds are rather 
close-by and yet they already have negative line-of-sight veloci- 
ties. For example, the range of GSR velocities at x' — —5 kpc 
ranges from -20 to20kms~^, while at —4 kpc the range is -50 to 
— 20knis~^. It is somewhat unlikely that these clouds would be 
able to wrap around to the other side of the Galactic centre with 
such infalling velocities. The near method in Section IZTI applied 
to Complex C gave similar but slightly smaller distances. Both the 
current treatment and the near method have likely given distances 
which are too small to be plausible. We deduce that Complex C is 
not a stream lying in a plane through the Galactic Centre. 



3 CORRECTION OF RADIAL VELOCITIES 



The method discussed in Section [Z2l was revisited, but with a cor- 
rection now applied to the velocities in order to account for our 
offset from the Galactic centre. Whereas previously the heliocen- 
tric line-of-sight velocities corrected to the GSR were used in the 
energy equation, we now corrected these iteratively to actual Galac- 
tocentric radial velocities. We first calculated, using line-of-sight 
velocities, the zeroth order distances to the stream clouds. With 
the three-dimensional configuration of the stream then known to 
lowest order, we could transform to planar coordinates with the 
suppressed axis corresponding to the normal to the stream plane 
as defined by the three fiducial points (see Section [Z2l for defini- 
tion). We discounted the third dimension — vertical offsets from 
the stream plane — on the basis that these deviations were small 
relative to both components in the plane and overall distance of the 



stream from the Galactic centre. We then know the approximate 
direction in which a particular point in the stream is moving, as- 
suming no motion into or out of the plane. This assumption was 
reasonable given the size of the vertical offsets. The component of 
the true velocity projected along our line of sight gives the observed 
velocity vi, since radial velocities are taken to be positive if they are 
outwards, hence an appropriate de-projection of vi gives the 'true' 
velocity. This can then be re-projected along the radial vector from 
the Galactic centre to the stream point to give the corrected radial 
velocity. 

We found the line-of-sight velocities to have been a slight 
under-estimate on the Galactocentric radial velocities for the Mag- 
ellanic Stream. The maximum change in Galactocentric distance as 
a result of this velocity correction is approximately 3.2 kpc for the 
farthest point, and the smaller distances are less affected, on aver- 
age by less than 1 kpc. We conclude that the size of the corrections, 
which are presented in Table|2l are small in comparison with errors 
that might result from the choice of method to derive the distances. 

By calculating the 'true' velocity vector, we were also able 
to calculate the angular momentum vector at each point along the 
stream. In Figure [4] the magnitude of the a::/-component of the 
angular momentum is plotted as a function of distance along the 
stream, where the xf axis is parallel to the normal to the stream 
plane, and nearly parallel to the Galactic x axis. 

In order to determine the component of velocity along the 
stream, we require the angle between the stream and the vector 
transverse to the radial vector. This can be determined from Figure 
[3] Since the angular momentum depends on the cotangent of this 
angle, the errors on the angular momenta are derived from assum- 
ing a ±2° error in this angle. Near the perigalacticon, the angle be- 
comes close to 90° . For such points, the associated errors become 
very large; these points have not been plotted. Figure |4] shows that 
the angular momentum along the stream is consistent with being 
constant within the errors. 



4 EFFECTS OF THE GALACTIC POTENTIAL 

We now address the issue of the choice of potential used. Until now, 
we have assumed a simple logarithmic potential giving a constant 
circular velocity. It is instructive to check whether this is an over- 
simplification and we should instead be using a potential which has 
a power-law dependence on the Galactocentric distance. As with 
the radial velocity correction above, the baseline for our check is 
the method given in Section [2J2] The following equation results 
from using properties of the H point in combination with a rear- 
rangement of equation ^ with tp{r) oc r~": 



1 + 



vl 



^1 -X-^-^il-X' 

Ho 



= 



(17) 



where X = r/rn and all other variables have the same definitions 
as before. 

When the constant a is set to a small positive quantity, a small 
negative correction results for the Galactocentric distances. For ex- 
ample, a — 0.01 shifts the distances found with the logarithmic 
potential down by approximately 0.7 kpc and a = 0.1 shifts the 
distances by approximately 7 kpc. 

It is instructive to find the range of values of a that is consis- 
tent with acceptable values of the Galaxy's mass out to 100 kpc, 
A/ioo- Imposing the total mass within this radius and a circular ve- 
locity at the Sun automatically imposes the steepness of the density 
profile. Since Tp{r) oc r~" and —r^dTp/dr = GM, it follows that 
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GM(r) — v^R"r^ ". If we ask, for example, that a be consistent 
with 0.8 X 10^^ Mq < Mioo < 2.0 x 10^^ Mp which is consis- 
tent with values given bv lKulessa & Lvnden-Belll ( Il992h . then a is 
restricted to values between 0.14 and —0.24. Employing a smaller 
mass for the Galaxy would raise the value of a. If, for example, the 
mass range given above were applied at 200 kpc, then the value of 
corresponding a would range between 0.33 and 0.04. The actual 
range is narrower when physical demands are incorporated, such 
as that the Magellanic Clouds should not be at ridiculously small 
distances and that we should perhaps expect the circular velocity 
to taper off gradually between the solar position and the distance at 
which the farthest point on the stream resides. We would, for exam- 
ple, expect a to be no greater than 0. 1 in order that the distance to 
the Magellanic Clouds are not too discrepant with observationally 
determined values. 



5 RELAXING THE ANGULAR MOMENTUM 
CONSTRAINT 

In previous sections, we have assumed all of the material compris- 
ing the stream to be following the same orbit, with the same value 
of angular momentum and energy being applicable along the entire 
length of the stream. We now investigate the effect of relaxing this 
constraint. If the material forming the Magellanic Stream had been 
stripped from the Clouds at the last passage through perigalacticon, 
the material must have undergone a change in energy and angu- 
lar momentum in order for the observed streaming shape to have 
formed. There is an observed lag of approximately 95° of the end 
of the stream relative to the position of the Magellanic Clouds. 

In the course of a tidal tearing, the outer parts of the victim are 
torn off by the strong tides close to pericentre. For a small victim, 
the tidal debris will have almost the same pericentre as the remnant 
of the victim. Thus the energies and angular momenta of the tidal 
debris are related by the condition of all having the same pericentre 
distance. An alternative interesting way of seeing this is to imagine 
that near pericentre, where the tearing occurs, it is sufficient to treat 
the potential under which the debris moves as that of two bodies 
rotating about one another with angular velocity Qp . In reality the 
rotation rate varies, but it will come to a maximum at pericentre so 
it will not be a great error to treat it as constant near there. Now 
an orbit in the field of two bodies rotating around one another at 
constant angular velocity Up will have a constant Jacobi integral 
(or energy relative to the rotating axes). The quantity e — flp^ will 
be conserved as debris leave the victim. This implies that those that 
gain most angular momentum Ah from lagging just behind the vic- 
tim will gain energy rip Ah. This argument leads to essentially the 
same conclusion as the demand that the perigalactica of all the de- 
bris be the same. Hence when we deal with debris spread along the 
orbit, we change both energy and angular momentum in proportion 
and use Ah as our parameter. 



5.1 Method 

We continue to use a generalised potential of the following form in 
the energy equation, now with an explicit constant: 

V'M = J- (18) 

Differentiating the potential and using the fact that the square of the 
circular speed at radius r is given by —rdtp/dr gives the equation 
lie (r) — aK/r". Since we have both a value for the circular speed 



at the location of the Sun and its distance from the Galactic centre, 
K can be calculated given a value for a; the potential at a given ra- 
dius is then known. We can rearrange the energy equation by using 
u = 1/r to obtain: 



u' +u = T^ (V'(w) + e) 



h^ 



(19) 



where the prime denotes differentiation with respect to the az- 
imuthal angle <jf>. Using equations l[9j, l |18l l and l |19l l for the peri- 
and apogalacticon gives the following equation for the angular mo- 
mentum h: 



2^(MpcH 



(20) 



Keeping the assumption that the Clouds and stream lie in a planar 
configuration, the orbit of the Cloudy can be computed for a given 
perigalacticon and apogalacticon. The constant a in the potential 
was set to 0.01. Note that it is not strictly necessary to have an 
absolute value of K if we only require relative values of angular 
momenta. 

The true current position of the Clouds was estimated to be 
approximately 0.7 rad from the perigalacticon that the Clouds have 
just passed through; we refer to this angle as the 'offset' angle. Note 
that more than half of the stream material lies on the other side of 
this perigalacticon. In the orbital computation, we know the time 
it has taken for the Clouds to cover the distance from the initial 
perigalacticon — where we assume the material for the stream was 
originally stripped — to their current position. During the same 
period in time, the end of the stream must have covered a smaller 
angle in the sky in order to have left the observable lag-angle of 
approximately 95°, despite having started at the same location. 



5.2 Results 

We took a fixed perigalacticon distance for the Clouds and investi- 
gated the effect that changing the distance to apogalacticon has on 
the total change in angular momentum that occurs along the stream. 
For rapo of between 80 and 180 kpc and a fixed rpori of 45 kpc, we 
find a simple relationship between r-apo and Ah/h through a fit of 
the following form: 



Ah 
IT 



19.4 kpc; 6~ 1.55, 



(21) 



where Ah now specifically refers to the absolute change in angular 
momentum from the Clouds to a point 95° downstream. The largest 
difference in angular momentum arises from the least eccentric or- 
bits. For an assumed apogalacticon of 120 kpc for the Magellanic 
Clouds, we find that Ah/h — 0.060, implying a 6% difference in 
the angular momentum between the extreme ends of the stream. 
For these choices of parameters, the Clouds would be positioned 
at r = 51.7 kpc and the end of the stream at 59.3 kpc. Changing 
only the perigalacticon to 44 kpc brings the Clouds slightly further 
in to 50.6 kpc and a slightly smaller value of Ah/h. We find that 
the distance in an unperturbed orbit to the same azimuthal position 
relative to the Magellanic Clouds as the end point is only of order 
0.5 kpc less than that for the perturbed orbit. 



•^ For the puiposes of this calculation, the Large and Small Magellanic 
Clouds — collectively referred to as 'the Clouds' — were taken as a com- 
bined single point at their approximate centre of mass. 
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5.3 Consistency Checks 

By allowing for a gradient in the angular momentum along the 
length of the stream, the angular position of the current perigalacti- 
con will be different between the Clouds and the end of the stream, 
even though, by design, the distance of closest approach to the 
Galactic centre will remain the same. When setting an offset angle 
for the Clouds, which is taken relative to the current perigalacti- 
con, it is important to remember that the location of perigalacticon 
gradually shifts along the stream. We have checked that its differ- 
ence for the extreme ends of the stream is merely 1 — 2° and hence 
there is little error in measuring the offset angle from the perigalac- 
ticon of the computed Magellanic Clouds' orbit, as opposed to an 
averaged value. 

Using the first set of parameters described in the section 
above, we have also checked the effect of altering the offset angle 
by a few percent to ensure that an error in determining this value 
from observational data and calculations described in Section |Z2l 
would not lead to a significant error in Ah/h. We find that the per- 
centage difference in this quantity resulting from a +10° change 
in the offset angle is tiny, on the order of a few hundredths of a 
percent. An under-estimation of the offset angle by 10° produces 
a greater error, on the order of a few tenths of a percent of Ah/h. 
This is as expected, given that changes in velocity over the same 
angle occurs more rapidly when closer to perigalacticon than when 
further away. The change in angular momentum that would be re- 
quired to keep the end of the stream far enough away therefore be- 
comes slightly larger. Even for the under-estimation case, however, 
the errors are still relatively small. 

One other point to note is that the concept of an H point, de- 
fined previously as the location on the orbital stream with the high- 
est radial velocity, is now void. This is because the stream material 
is not following the same orbit as that of the Clouds. A check to see 
whether the projected line-of-sight velocity at the tip of the stream, 
earlier coined the H point, is consistent with the observed value 
at that location served as an independent test of the validity of the 
approach used to exploit the lag-angle in calculating the angular 
momentum deviation. 



6 SUMMARY AND CONCLUSIONS 

We have applied a variety of geometrodynamical methods to 
streams of high- velocity clouds in an attempt to constrain their dis- 
tances. Some of these methods, including the near method, have the 
nice feature that distance determinations would be possible with- 
out the parameter dependency inherent in any simulations. We find, 
however, that the near method in particular becomes rather inaccu- 
rate when the angle subtended at the stream by Ro becomes small. 
At the expense of becoming more model dependent by assuming 
a halo form for the Galactic potential, we have improved on this 
accuracy. We then find distances to the entire length of the Magel- 
lanic Stream by assuming that the gas comprising the stream nearly 
follows a planar orbit about the Galactic centre. We find the peri- 
galacticon to be at 45 kpc and the maximum distance at the trailing 
end to be 75 kpc. The heliocentric distances are not very differ- 
ent from the Galactocentric distances, varying typically by 1 kpc 
or less. The centre-of-mass of the Magellanic Clouds would be at 
47 ± 5 kpc from us, according to this method; this is consistent 
with various obser vationa l ly det ermined values for the Large Mag- 
ellanic Cloud (see I Alvesl J2004r) for a review). Application of our 
methods to Complex A and C did not yield very promising results 



under the assumptions that were made. This indicates that our hy- 
pothesis, namely that the streams were nearly following an orbit 
around the Galactic centre in a well-defined plane, was not valid 
for these two streams of high- velocity clouds. 

In this paper, we have tested the idea that streams of HVCs 
follow orbits. However the assumption that the orbit lies in a plane 
that passes through the Galactic centre is rather restrictive when 
applied to nearby streams that do not pass close to the Galactic 
Pole. It is possible to remove the restriction of a planar orbit, pro- 
vided a model of the Galaxy's potential is given. We find that the 
transverse velocity of the stream can be deduced from the run of 
radial velocities as a function of angle along the stream, provided 
the distance to one point of the stream is assumed. Orbits can then 
be computed in the given potential. Both directions in the sky and 
radial velocities are then computed and these are compared with 
the observations. Such a method and an example of its applica- 
tion to Co mplex A and the stellar Orphan stream has already been 
presented iJin & Lvnden-Bel]||2007h . We note that such methods 
of deriving the tr ansverse velocity comp onent were used in a re- 
lated problem by iFeitzinger et alj ( Il977h . but probably date back 
to an earlier century. We are also exploring the possible existence 
of stream associations. Following on from these works, we hope, 
eventually, to give a definitive answer to the question of whether 
high- velocity cloud streams follow orbits around the Galaxy. 
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dle: line-of-sight velocity along the orbit. Bottom: the angle i?^ subtended 
at the stream by Ro projected onto n. Initial conditions: Galactocentric 
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V = (—151, —143, 0.2) km s~^. The time at which the Z point is reached 
is denoted by the dot-dash line. 
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Figure A2. The variables X and Y of equation (8) have been plotted against 
one another for the orbit shown in figure IaT] The straight solid lines show 
the results from an exact orbit (where, for clarity, individual points are not 
shown). The crosses give the results for the same orbit, but where the direc- 
tion cosines have had a random scatter (of maximum 0.01) applied, before 
the orbit is smoothed by the code. The few points along the lower solid line 
indicate the behaviour of the X and Y variables for positions near the Z 
point. 



clearly visible in Figure lAZI both when using the exact orbital data 
and when some scatter has been applied to the direction cosines. 
The distance to the Z point, dz, can be calculated from the noisier 
plot to within a few per cent of the correct value. 



APPENDIX A: EXAMPLE OF THE NEAR METHOD 

We have tested the near method (see Section lOT ) on simulated or- 
bits of varying orientations and distances in the logarithmic poten- 
tial. The constraint on &d given in Section [ZT] was deduced from 
these test cases; one working example is shown here. The normal 
to the orbital plane is given by n = (-0.687, 0.725, 0.047). The 
orbit is nearly polar but is otherwise inclined at approximately 43° 
to the Galactic x axis. 

Figure lATI shows how the heliocentric distance to the orbit, the 
line-of-sight velocity and the angle &d = arcsin(Ro.n/d) change 
as a function of time. The orbit shown has a length of 130°. The 
singular behaviour of the variables X and Y near the Z point is 



